Theory of magnetic structure in layered iridates: spin-orbit band or Mott insulators 
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In iridates, owing to strong spin-orbit coupling, the effective total angular momentum j = 1/2 
band is well separated from the rest of the bands. In particular S^IrCU (Sr-214) and Sr3lr207 
(Sr-327) are magnetic insulators formed in j — 1 /2 bands leading to an interesting proposal named 
spin-orbit Mott insulators. However, given that there is an even number of Ir atoms per unit cell due 
to a staggered rotation of octahedra and even number of layers, the non-interacting system could 
be a spin-orbit band insulator, questioning the Mott character in these materials. To answer the 
question on the nature of the insulating phases, and different states in various layered structures, we 
study a Hubbard model with a tight-binding spectrum designed for each Sr-214, Sr-327 and SrlrOs 
(Sr-113). A canted antiferromagnet (AF) is found in Sr-214 which is deep in the insulating phase 
close to strong coupling limit, while a collinear AF with c-axis moments is realized in Sr-327 near 
the spin-orbit band insulator. In contrast, Sr-113 is semimetallic. The origin of such a disimilarity is 
explained, and implications of our results in relation to possible high temperature superconductors 
in doped iridates are further discussed. 

PACS numbers: 71.30.+h,71.70.Ej,75.30.Kz 



Introduction - Spin-orbit coupling (SOC) has often 
been either ignored or treated perturbatively in corre- 
lated electronic systems such as 3d- and 4d-orbital mate- 
rials. However, with heavier atom systems, it becomes an 
important ingredient in determining the ground state of 
these materials and leads to novel phases such as topo- 
logical insulators and spin liquids. In particular, when 
the SOC strength is comparable to the electronic inter- 
actions, understanding their interplay becomes challeng- 
ing. Ir-oxides (iridates) with 5d-orbitals offer such a play- 
ground to investigate their combined effects. 

It has been suggested that the effect of Hubbard inter- 
action in iridates is amplified due to strong SOC via the 
narrowing of bandwidth, leading to an insulating state in 
some layered perovskitefTH2rj] and pyrochlore[271 [55] iri- 
dates. This state, called a spin-orbit Mott insulator, was 
first reported in single layered perovskite Sr2lr04 (Sr- 
214). [6l [3 |9] An anisotropic Heisenberg spin model was 
derived at large U limit [5J, and it was further proposed 
that anisotropic terms can be gauged away leading to an 
isotropic Heisenberg model of j = 1/2, reminiscent of the 
spin-1/2 undoped high temperature cuprates. Using this 
simple mapping, the idea of d-wave high temperature su- 
perconductivity in doped iridates was proposed. |14j 

In this paper, we first question if the iridates are in- 
deed Mott insulators and further investigate whether the 
mapping to the isotropic Heisenberg model for j — 1/2 is 
proper. To address these questions, we study a Hubbard 
model taking into account a tight binding Hamiltonian 
for the full unit cell. One important fact to notice is 
that in Sr-214 and Sr-327, there is an even number of 
Ir atoms due to the staggered IrOe octahedra rotation 
and four layers per unit cell. Thus according to the band 
theory, the non-interacting system could be a spin-orbit 
band insulator. 
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FIG. 1. [Color online] Phase diagram for the magnetic phases 
of the single layer, bilayer and 3D systems with the accom- 
panying metal-insulator transition. The critical value of U 
of the magnetic transition for each system is U c /t = 1.07, 
1.72 and 1.98 respectively. For the Sr-214 and Sr-113 system, 
since the underlying band structure is metallic, the transition 
turns the metal into a magnetic metal (MM) before it devel- 
ops a gap as the ordering increases. The ordered phases in the 
Sr-113 system is a non-coplanar canted AF (N-CAF), where 
the ferromagnetic moment points along the crystal c-axis. In 
a scenario where the Hubbard interaction is 1.73t ~ 0.5eV 
(setting t ~ 300meV), the stars represents where each sys- 
tem should be on the phase diagram. 

We found that Sr-327 is indeed a band insulator at 
small U , which becomes a magnetic insulator as U in- 
creases. In contrast, Sr-214 is a metal at small U, which 
then turns into a magnetic metal and a magnetic insu- 
lator subsequently, while Sr-113 remains semimetal even 
at relatively large U as shown in Fig. 1. 

A collinear AF with moments oriented along the crys- 
tal c-axis is found in Sr-327, while a coplanar canted AF 
is realized in Sr-214 as shown in Fig. 2. While the mag- 
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FIG. 2. [color online] Theoretical results of magnetic order- 
ing in Sr-214 (a) and Sr-327 (b) in the insulating phase. Each 
layer has two different atoms because of the staggered rota- 
tion of the octahedra which we label blue and red. (a) The 
coplanar canted AF state is shown with the spin configura- 
tion displayed for each layers. The black arrow represents 
the direction of the ferromagnetic moment within the plane, 
showing the up-up-down-down structure, (b) The favored 
collinear AF state with moments pointing along the c-axis is 
shown where the second set of bilayer have its spins flipped 
from the first set. 



netic ordering is certainly an effect of the interaction, it 
is not sufficient to conclude that iridates are Mott insu- 
lators in the strong U limit, where the charge gap due 
to electronic interaction should remain finite even when 
the magnetic order disappears. Can we distinguish weak 
U vs. strong U by studying their magnetic ordering pat- 
terns, which then act as a guide to understanding Mott 
insulators? We begin this study with a tight binding 
model. 

Tight-binding model - A consistent description of the 
Ruddlcsdcn-Popper series of strontium iridates can be 
achieved using a tight-binding model with on-site Hub- 
bard interaction. Nearest in-plane and intra-bilayer hop- 
pings (represented by t p (t' p ) and t z (t' z ) between same 
(different) j z state, respectively), next-nearest (t n ), next- 
next-nearest (t nn ) neighbor hoppings, and nearest-inter- 
(bi)laycr hoppings {t u t id , t lz , t iy and t i x ) in the j — 1/2 
basis are included. Since both the Sr-214 and Sr-327 
structures have four distinct layers per unit cell, and each 
layer has two different Ir atoms due to a staggered Ir06 
rotation, it is convenient to introduce the following spinor 
basis: V£ = (c£ t B , cgf, cgf, cgf), where a = 1,...,4 is 
the layer index, B(R) is the sublattice label representing 
alternating IrOg rotation, and arrows denote j z = ±1/2. 

Since we only consider hoppings between neighboring 



layers, the Hamiltonian can be written in general as H = 
Efc ^f H f^k + UJ2i 4f c 4 c 4 c iT> where 

/ A k B k e ifc * c C|\ 

TraP _ B k A k C k m 

Hk ~ C\ A k B k ' (1) 

\e~^ c C k B\ A k J 

and U is the on-site Hubbard interaction in j-basis. Note 
that in addition to bandwidth effects, the Hubbard in- 
teraction U is renormalized by a factor of 1/3 from U 
of the orbital basis, due to basis change from the orbital 
basis to the j basis (i.e., U = U/3). Here, A k , B k , and 
C k are 4x4 matrices with A k containing intra-layer hop- 
pings only, B k contains hoppings between layers 1 and 

2 (3 and 4) and C k contains hoppings between layers 2 
and 3 (4 and 1 where c is a c-axis lattice contant). Since 
A k is within a given layer, it is the same for both the 
single-layer and bilayer systems. 

The main difference between the single-layer and the 
bilayer system is the coupling between layers 1 and 2 (and 

3 and 4), Bk- In the single-layered version, this coupling 
is weak and each atom is coupled with 4 atoms of the next 
layer (because of the c-axis stacking) whereas it is strong 
in the bilayer Sr-327. Note that R (B) is sitting on top 
of B (R), which differs from Sr-214. This is responsible 
for the different magnetic ordering patterns for these two 
layered iridates, which we will discuss further below. 

For both the single layer and the bilayer system, the 
coupling between layers 2 and 3 (and 4 and 1) are nearly 
identical due to the similarities between the inter-layer 
structures with layer 3 being displaced from layer 2 by 
half a planar lattice vector (110 or 110). We therefore 
allow the Ck matrices for both system to be identical. In 
Sr-214, the small difference between Bk and Ck is due to 
the difference in displacement between neighboring lay- 
ers. Between layers 1 and 2, the displacement is half the 
lattice vector (110) whereas it is half the lattice vector 
(110) between layers 2 and 3. All the forms of Ak Bk 
and Ck which include the dispersions and hopping pa- 
rameters derived using Slater-Koster theory are given in 
the supplementary materials. 

Mean-Field Theory - Given the tight binding Hamilto- 
nian, we then treat the Hubbard interaction at the mean- 
field (MF) level to find possible magnetic orderings. Per- 
forming a MF decoupling in the magnetic channel, the 
interaction is written as, 

Him = ~{SL + Sf y + Si) + I ^(ntt + n u ) (2) 

i 

where we take Sf a = 2mi a Si a — mf a and compute m,i a 
self-consistently. There are 24 order parameters from 4 
layers, 2 sublattices, and 3 directions of m for Sr-214 
and Sr-327. In Srlr03 (Sr-113), there are 12 order pa- 
rameters for 4 sublattices and 3 directions of m. As U 



3 



increases, there is a second order phase transition from 
metal/semimetal/band insulator into a magnetically or- 
dered phase as shown in Fig. [I] Below we discuss some 
details of phase diagram for each system. 

For Sr-214, the magnetic ordering occurs at U c = 1.07t, 
where t is the magnitude of the Slater-Koster hopping 
tdd-n and used as the energy scale throughout the remain- 
der of the paper. The canted AF is the ground state. In 
addition, the direction of the net ferromagnetic compo- 
nent is also determined to be either in the (110) direc- 
tion or the (110) direction in an up-up-down-down (110) 
or up-down-down-up (110) pattern along the c-axis as 
shown in Fig. [2] This magnetic ordering pattern re- 
mains unchanged in the large [/-limit, and is consistent 
with the experimental results obtained by resonant x-ray 
scattering techniques. [TJ [31 [SI IH] 

For Sr-327, when U < U c = 1.72t, one has a spin-orbit 
band insulator. It then becomes a magnetic insulator, 
where the charge gap increases gradually across the sec- 
ond order transition. Note that there are five B Ir-atoms 
surrounding R Ir-atom in Sr-327. The competition be- 
tween canted AF and collinear AF strongly depends on 
the in-plane and z-axis bilayer hoppings. We found that 
only when the ratio of t' z /t z is close to the ratio of t' p /t p , 
the degeneracy of canted AF and collinear AF persists. 
We estimate t' z /t z w 3 t' p /t p for Sr-327 based on Slater- 
Koster theory, leading to a collinear AF state along the 
c-axis favored by the bilayer coupling. Such a state has 
been recently confirmed in experiments [2"U1 |2"T1 [231 l2"Sj . 
The canted AF state can become a lower energy state 
when the inter-layer couplings (Cj.) are enhanced, but 
only in an unrealistic limit, or when the ratio between 
t' z /t z and fp/tp is very close to 1. Due to the two sets 
of bilayers per unit cell, 2 different configurations with 
almost degenerate energies are possible in the collinear 
AF ground state. One (with a slightly lower energy) is 
shown in Fig. [2] (b), and another spin configuration is 
equivalent to having the opposite spins on the 2nd bi- 
layer set. These two configurations were experimentally 
observed in ref. |23j . 

The difference in energy between the collinear AF state 
and the coplanar canted AF state can be used here to es- 
timate the spin gap reported by resonant inelastic x-ray 
scattering (RIXS) techniques. For the set of parameters 
given in the supplementary material, the difference in en- 
ergy between the 2 states per Ir atom for U/t — 2 is of 
the order of 3 x 10 _4 i. Approximating the Slater-Koster 
hopping \tddTrl — 300meV, we get a spin gap which is ap- 
proximately 0.1 meV, two orders of magnitude smaller 
than the 90meV gap[21 . This result seems to indicate a 
small spin gap, which may be hidden within the resolu- 
tion of the RIXS data. 

The same mean-field approach can be used on the or- 
thorhombic perovskite Sr-113 model described in [18]. In 
this case, because of the tilting of the octahedra (in addi- 
tion to the rotation) , a non-coplanar canted AF magnetic 



ordering occurs, but with the canted moment pointing 
along the z-direction instead of in the plane. The mean- 
field transition occurs at U c — 1.98£ in this case, and the 
Hubbard repulsion in iridates should be ~ 1.73£ corre- 
sponding to ~ 0.5eV. 

Large-U limit spin model - We found that the mag- 
netic structures do not change with a large U in our 
MF study. To understand how the above weak cou- 
pling results are related to the strong coupling ap- 
proach, we derive a spin model in the large U limit. 
Within a second order perturbation theory, the Heisen- 
berg, Dzyaloshinskii-Moriya (DM) and anisotropic ex- 
change terms are obtained as follows. For a general time 
reversal symmetric hopping between two sites given by 
(t a p + ■ Vap)4,a0 C j,P<r' > the spin-model is 

<ij>aP 

(3) 

Here J?/ = £(t^ - v^), Dff = §(2t a pv a0 ), and the 

anisotropic exchange Tf^ ab = ^(2u*^u^) where a,b are 
spin components and a/3 are layer indices. Note that 
these exchange terms depend not only on site i, j but also 
on layer index and spin component. The in-plane nearest- 

4(t 2 -t' 2 ) 

neighbor exchange can be written as J BR = — v - u p , 

Der = —jf-^-z and zz = and are the same for 
Sr-214 and Sr-327, for a = 1, 4 leading to a degeneracy 
between canted AF and collinear AF if the inter-layer 
terms are ignored. 

For Sr-214, the exchange terms between layers 1 and 2 
is found to be J 12 B = ^, J 12 R = ±(t 2 d - t 2 iz - t% -t 2 x ), 

^ BR ~ JJ (-t-tidtix > ^-tidtiyitidtiz) with 

/ t 2 t t ±t- t \ 

r 12 - — t t t 2 +t- t- 

\±t t ±t- t t 2 I 

where U represents inter-layer hopping integral between 
B and B (or R and R), while ti d between B and R. The 
spin cr Q -dependent hopping term between B and R are 
denoted by U a where a = x,y, z. The ± denotes the 
sign change in the two a x - and <x y -dependent hoppings 
between B(R) in layer 1 and the two R(B) in layer 2 it 
connects to. Because of the sign changes, these terms 
have effectively no effect if the magnetic unit cell is the 
same as the lattice unit cell. The exchange terms between 
layers 2 and 3 are identical except for the y-component of 
the DM-term where D BR — —D B ^ and the xy- and yz- 
components of the anisotropic exchange where r^ 3 j^ 1 ' = 
-i BR ana i BR — i BR . 

The isotropic term by itself would favor a collinear AF 
state, but it is largely frustrated due to the fact that each 
atom is coupled with 4 atoms of each neighboring layers. 
The same applies to the DM-vectors but the anisotropic 
exchange term T xy on the other hand favors the canted 
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AF state in an up-up-down-down configuration. Due to 
the reduced effect of the other terms because of the lat- 
tice structure, the anisotropic exchange is large enough 
to have the canted AF be the ground state. It has been 
shown 8] that the Hund's coupling and tetragonal dis- 
tortion also favor the canted AF state, further stabilizing 
this state. 

For Sr-327, the in-plane spin exchange are identical to 
those for the single layer. In addition, the bilayer ex- 
change are given by J B 2 R = ±(t 2 z - t' 2 ), D B 2 R = e^t^z 
and T l RRzz = jjt' 2 . Note that unless t' z /t z — t' p /t p , the 
DM interaction in-plane and out-of-plane are frustrated, 
i.e. it is impossible to configure the spins such that the 
optimal canting is achieved both in-plane and within the 
two layers. Thus, the collinear AF state is the ground 
state, due to the opposite rotation of the octahedra be- 
tween adjacent layers, consistent with the bilayer spin 
model results discussed in Ref. [20] . 

Adding the inter-layer hoppings, the c-axis stacking fa- 
vors an AF pattern between the same atom of layers 2 
and 3 as seen in Fig. [2] but the ferromagnetic pattern is 
very close in energy, leading to the two possible arrange- 
ments discussed previously. 

Discussion and summary - We offer a systematic study 
for different layered iridates using a Hubbard model that 
captures a rich phase diagram; metal to magnetic metal 
to magnetic insulator for Sr-214, band insulator to mag- 
netic insulator for Sr-327, and semimetal to magnetic 
semimetal to magnetic insulator for Sr-113, as U in- 
creases. The ground state magnetic ordering patterns 
are sensitive to the lattice structure. In Sr-214, the inter- 
layer isotropic Heisenberg and DM interactions are frus- 
trated, because Ir atom in adjacent layers have zero ef- 
fective field due to their position. This effect, combined 
with the anisotropic exchange term, selects the "up-up- 
down-down" pattern for the net ferromagnetic moment 
pointing in (110)- or (llO)-direction. On the other hand, 
the bilayer lattice structure frustrates the DM interaction 
of the 5 neighbors of each Ir atom making the collinear 
AF state along the c-axis the preferred configuration in 
Sr-327, in contrast with Sr-214. However, the canted AF 
is close in energy, making the spin dynamics, tempera- 
ture and magnetic field dependence interesting subjects 
for future study. 

For all the iridates studied here, magnetic ordering 
structures are identical to those for large U limit ob- 
tained by the spin model, implying that their Mott char- 
acter cannot be corroborated by the magnetic ordering. 
However, Sr-214 is far from its metal-insulator transi- 
tion, deep in the insulating phase, while Sr-327 is nearby 
a small gap band insulator at U ~ Q.beV suggesting that 
the magnetic ordering does not support a Mott insula- 
tor. We suggest that ARPES study below and above the 
transition temperature to test our proposal. 

A comparison to cuprates deserves some discussion. 
Given that j = 1/2 spin model resembles the high tem- 



perature cuprates of spin 1/2, a possibility of achieving 
d-wave high temperature superconductivity in doped iri- 
dates became appealing. However, the essential differ- 
ence is that the j = 1/2 wavefunction is made of equal 
mixture of the ti g orbitals, leading to an isotropic s- wave- 
like wavefunction while the two dimensional d x 2_ y 2 or- 
bital is a basic ingredient of the cuprates. This results 
in a strong bilayer coupling in Sr-327 at the origin of the 
spin-orbit band insulator at small U regime, making this 
system distinct from bilayer cuprates where the bilayer 
coupling can be treated perturbatively. Along similar 
lines, the anisotropic exchange terms are not negligible 
and the proposed d-wave pairing in doped iridates will be 
unlikely achieved, even though superconductivity is still 
an open possibility. 
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SUPPLEMENTARY MATERIAL 
Hopping parameters and dispersions 

The intra-layer hopping matrix is given as 

A k = 4 + efr x + e' k ad r y a z , (4) 

where r represents the sublattice degree of freedom and 
a represents the spin. The dispersions obtained using the 
symmetry of the hopping parameters are given as 

£fc = 4t„ cos(k x ) cos(fcy) + 2t nn (cos(2k x ) + cos(2k y )) 
e% d = 2t p (cos(k x ) + cos(ky)) 
e' k ad = 2t' p (cos(k x ) + cos(k y )), 

with the hopping parameters listed in Table I. 
The first inter-layer hopping matrix given by 

Bik = e b k + e b k d T, x + e b k z r y a z + e b k y T y a y + e b k x T y a x , 

B 2k = t Z T x + t' z T y (T z (5) 

where B\ k is for the Sr-214 system with weak inter-layer 
hoppings and B 2k is for the Sr-327 system with strong 
bilayer coupling. The dispersions for the inter-layer cou- 
plings in matrix B\ k are given by 



4 


= 2ti cos((fc x 4- 


-ky)/2) 


hd 

e k 


= 2t id cos((k x 


~ ky)/2) 


bz 
€ k 


= 2t iz cos((fc x 


~ ky)/2) 


bu 
C k 


= i2t iy sm((k x 


- ky)/2) 


bx 
c k 


= i2t ix sin((fc K 


~k v )/2), 



whereas B 2k doesn't have a dispersion since the B(R) 
atom from layer 2 sits on top of the R(B) atom of layer 
1. 

Finally, the second inter-(bi)layer hopping matrix is 
given by 

C k =4 + f*T w + e c k z T y a z + e c k y T y a y + e c k x T y a x , (6) 

with the dispersions given by 

e c k = 2Ucos((k x - k y )/2) 
ef = 2t ld cos{(k x + k y )/2) 
ef = 2t iz cos((k x + k y )/2) 
e c k v = -i2t iy sm((k x + k y )/2) 
e c k = i2i lx sin((fc a; + k y )/2). 

The hopping parameters for B\ k and C k are listed in 
Table I. 

All the hopping parameters were derived using Slater- 
Koster hoppings for d-electrons with (tddir,tdd<T,tdds) = 
( — 1, 3/2, 1/4) for nearest neighbor hoppings, and setting 



TABLE I. List of intra-layer hopping parameters and list of 
bilayer and inter- (bi) layer hopping parameters in unit of t 



tp 




t' 

V 










-0.57 




-0.10 




0.005 




-0.01 


t z 


a. 


ti 




t / z 


tiy 


t i x 


-0.46 


-0.23 


-0.029 


-0.0275 


-0.0135 


-0.0095 


0.0095 



the ratios between next-nearest-neighbor and nearest- 
neighbor hoppings to be r = 0.2 and between next-next- 
nearest neighbor and nearest-neighbor to be r' = 0.02. 
We also used the ratio between inter-(bi)layer hopping 
to nearest-neighbor hopping to be r" = 0.1. 

Note that bilayer hoppings t z and t' z are relatively sim- 
ilar to the in-plane hoppings of t p and t', respectively. 
As we discussed in the main text, this is inherited from 
the character of j = 1/2 wavefunction containing equal 
weights of d xz , d yz in addition to d xy orbitals. This com- 
bination of j = 1/2 wavefunction also leads to an ex- 
tremely small next nearest neighbor hopping t n in the 
plane. There are three orbital hopping contributions 
along the diagonal directions between B to B (or R to 
R). One is from d xy to d xy orbitals whose hopping in- 
tegral is positive, and it is cr-type. On the other hand, 
d xz to d xz and d yz to d yz orbital hopping integrals are 
7r-type and comes negative. These two type hopping in- 
tegrals nearly cancel each other, and leave only a small 
strength of t n . The octahedra rotation makes the d xy 
hopping even smaller, which decreases t n further. 

Charge gap and magnetic moment 

Given the set of tight-binding parameters listed above, 
the Sr-327 system displays an underlying insulating band 
structure shown in Fig. 3 below. The direct gaps near X 
and M depend mainly on the in-plane and bilayer hop- 
pings (t p) t' p , t z and t' z ), and the overall dispersion also 
depends on the next-nearest in-plane neighbor hopping 

(tn)- 

While bilayer hoppings are necessary to open a direct 
gap, a large value of t n could push the bottom of the 
band near M further down, and push up the top of band 
near X further up in such a way that small electron and 
hole pockets at the X and M points could appear. These 
band touchings occurs when t n = 0.047 assuming that 
the other parameters are fixed. For the reason discussed 
the above, t n greater than 0.01 would be unrealistic and 
thus the insulating structure will be robust for the small 
U regime. 

We also note that the magnetic moment of Sr-327 is 
about 0.03/LtB@], which is close to our order parameter at 
U = 1.73i in Sr-327. This small order parameter should 
not affect the charge gap much, but the charge gap re- 
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FIG. 3. Underlying band structure for the Sr-327 model. This 
system is a "spin-orbit band insulator" at small U and a mag- 
netic transition occurs which opens the band gap further as 
U increases 

ported by a recent ARPES is of the order of 100 meV[2"5]. 
Therefore, this sizeable charge gap must originate mainly 
from the underlying band gap, especially because the cor- 
relation is not strong in Sr-327 as we shown in the current 
study. 



